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Abstract 

We present preliminary measurements of the CP-violating asymmetries in the decay /o(980)(— > 
The resuhs are obtained from a data sample of 209 x 10^ T(45) BB decays col- 
lected with the BABMl detector at the PEP-II asymmetric-energy B Factory at SLAC. From a 
time-dependent maximum-likelihood fit we measure the mixing-induced CP violation parameter 
S = -0.95lo;23 ± 0.10 and the direct CP violation parameter C = -0.24 ± 0.31 ± 0.15, where the 
first errors are statistical and the second systematic. 



Submitted to the 32 International Conference on High-Energy Physics, ICHEP 04, 
16 August— 22 August 2004, Beijing, China 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



The BABAR Collaboration, 

B. Aubert, R. Barate, D. Boutigny, F. Coudcrc, J.-M. Gaillard, A. Hicheur, Y. Karyotakis, J. P. Lees, 

V. Tisscrand, A. Zghichc 

Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, France 

A. Palano, A. Pompili 
Universitd di Bari, Dipartimento di Fisica and IN FN, 1-70126 Bari, Italy 

J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu 
Institute of High Energy Physics, Beijing 100039, China 

G. Eigen, I. Ofte, B. Stugu 
University of Bergen, Inst, of Physics, N-5007 Bergen, Norway 

G. S. Abrams, A. W. Borgland, A. B. Breon, D. N. Brown, J. Button-Shafer, R. N. Cahn, E. Charles, 

C. T. Day, M. S. Gill, A. V. Gritsan, Y. Groysman, R. G. Jacobsen, R. W. Kadel, J. Kadyk, L. T. Kerth, 
Yu. G. Kolomensky, G. Kukartsev, G. Lynch, L. M. Mir, P. J. Oddone, T. J. Orimoto, M. Pripstein, 

N. A. Roe, M. T. Ronan, V. G. Shelkov, W. A. Wenzel 

Lawrence Berkeley National Laboratory and University of California, Berkeley, CA 94720, USA 

M. Barrett, K. E. Ford, T. J. Harrison, A. J. Hart, C. M. Hawkes, S. E. Morgan, A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

M. Fritsch, K. Goetzen, T. Held, H. Koch, B. Lewandowski, M. Pehzaeus, M. Steinke 
Ruhr Universitdt Bochum, Institut fiir Experimentalphysik 1, D-44'if80 Bochum, Germany 

J. T. Boyd, N. Chevalier, W. N. Cottingham, M. P. Kelly T. E. Latham, F. F. Wilson 
University of Bristol, Bristol BS8 ITL, United Kingdom 

T. Cuhadar-Donszelmann, C. Hearty, N. S. Knecht, T. S. Mattison, J. A. McKenna, D. Thiessen 
University of British Columbia, Vancouver, BC, Canada VST IZl 

A. Khan, P. Kyberd, L. Teodorescu 
Brunei University, Uxbridge, Middlesex UBS 3PH, United Kingdom 

A. E. Blinov, V. E. Blinov, V. P. Druzhinin, V. B. Golubev, V. N. Ivanchenko, E. A. Kravchenko, 
A. P. Onuchin, S. L Serednyakov, Yu. L Skovpen, E. P. Solodov, A. N. Yushkov 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

D. Best, M. Bruinsma, M. Chao, L Eschrich, D. Kirkby, A. J. Lankford, M. Mandelkern, R. K. Mommsen, 

W. Roethel, D. P. Stoker 

University of California at Irvine, Irvine, CA 92697, USA 

C. Buchanan, B. L. Hartfiel 
University of California at Los Angeles, Los Angeles, CA 90024, USA 

S. D. Foulkes, J. W. Gary, B. C. Shen, K. Wang 
University of California at Riverside, Riverside, CA 92521, USA 



2 



D. del Re, H. K. Hadavand, E. J. Hill, D. B. MacFarlane, H. P. Paar, Sh. Rahatlou, V. Sharma 
University of California at San Diego, La Jolla, CA 92093, USA 

J. W. Berryhill, C. Campagnari, B. Dahmes, O. Long, A. Lu, M. A. Mazur, J. D. Richman, W. Verkerke 
University of California at Santa Barbara, Santa Barbara, CA 93106, USA 

T. W. Beck, A. M. Eisner, C. A. Heusch, J. Kroscberg, W. S. Lockman, G. Nesom, T. Schalk, 
B. A. Schumm, A. Seiden, P. Spradlin, D. C. Williams, M. G. Wilson 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, CA 95064, USA 

J. Albert, E. Chen, G. P. Dubois-Felsmann, A. Dvoretskii, D. G. Hitlin, I. Narsky, T. Piatenko, 

F. C. Porter, A. Ryd, A. Samuel, S. Yang 

California Institute of Technology, Pasadena, CA 91125, USA 

S. Jayatilleke, G. Mancinelli, B. T. Meadows, M. D. Sokoloff 
University of Cincinnati, Cincinnati, OH 45221, USA 

T. Abe, F. Blanc, P. Bloom, S. Chen, W. T. Ford, U. Nauenberg, A. Olivas, P. Rankin, J. G. Smith, 

J. Zhang, L. Zhang 

University of Colorado, Boulder, CO 80309, USA 

A. Chen, J. L. Harton, A. Soffer, W. H. Toki, R. J. Wilson, Q. Zeng 
Colorado State University, Fort Collins, CO 80523, USA 

D. Altenburg, T. Brandt, J. Brose, M. Dickopp, E. Feltresi, A. Hauke, H. M. Lacker, R. Miiller-Pfefferkorn, 
R. Nogowski, S. Otto, A. Petzold, J. Schubert, K. R. Schubert, R. Schwierz, B. Spaan, J. E. Sundermann 

Technische Universitdt Dresden, Institut fiir Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard, G. R. Bonneaud, F. Brochard, P. Grcnicr, S. Schrenk, Ch. Thicbaux, G. Vasileiadis, M. Verderi 

Ecole Poly technique, LLR, F-91128 Palaiseau, France 

D. J. Bard, P. J. Clark, D. Lavin, F. Muheim, S. Playfer, Y. Xie 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andreotti, V. Azzolini, D. Bettoni, C. Bozzi, R. Calabrese, G. Cibinetto, E. Luppi, M. Negrini, 

L. Piemontese, A. Sarti 

Universitd di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

E. TrcadwcU 

Florida A&M University, Tallahassee, FL 32307, USA 

F. AnuUi, R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, P. Patteri, L M. Peruzzi, 

M. Piccolo, A. Zallo 

Laboratori Nazionali di Frascati dellTNFN, 1-00044 Frascati, Italy 

A. Buzzo, R. Capra, R. Contri, G. Crosetti, M. Lo Vetere, M. Macri, M. R. Monge, S. Passaggio, 

C. Patrignani, E. Robutti, A. Santroni, S. Tosi 

Universitd di Genova, Dipartimento di Fisica and INFN, I-I6I46 Genova, Italy 

S. Bailey, G. Brandenburg, K. S. Chaisanguanthum, M. Morii, E. Won 
Harvard University, Cambridge, MA 02138, USA 



3 



R. S. Dubitzky, U. Langcncggcr 
Universitdt Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

W. Bhimji, D. A. Bowerman, P. D. Dauncey, U. Egede, J. R. Gaillard, G. W. Morton, J. A. Nash, 

M. B. Nikolich, G. R Taylor 

Imperial College London, London, SW7 2AZ, United Kingdom 

M. J. Charles, G. J. Grenier, U. Mahik 
University of Iowa, Iowa City, lA 522^2, USA 

J. Cochran, H. B. Crawley, J. Lamsa, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin, J. Yi 
Iowa State University, Am,es, I A 50011-3160, USA 

M. Biasini, R. Covarelli, M. Pioppi 
Universitd di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 

M. Davier, X. Giroux, G. Grosdidier, A. Hooker, S. Laplace, F. Le Diberdcr, V. Lcpeltier, A. M. Lutz, 
T. C. Petersen, S. Plaszczynski, M. H. Schune, L. Tantot, G. Wormser 

Laboratoire de lAccelerateur Lineaire, F-91898 Orsay, France 

C. H. Cheng, D. J. Langc, M. C. Simani, D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 

A. J. Bevan, C. A. Chavez, J. P. Coleman, I. J. Forster, J. R. Fry, E. Gabathuler, R. Garnet, 
D. E. Hutchcroft, R. J. Parry, D. J. Payne, R. J. Sloane, C. Touramanis 

University of Liverpool, Liverpool L69 72E, United Kingdom 

J. J. Back,^ C. M. Cormack, P. F. Harrison,^ F. Di Lodovico, G. B. Mohanty^ 
Queen Mary, University of London, El 4-NS, United Kingdom 

C. L. Brown, G. Cowan, R. L. Flack, H. U. Flaecher, M. G. Green, P. S. Jackson, T. R. McMahon, 

S. Ricciardi, F. Salvatore, M. A. Winter 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, 

United Kingdom 

D. Brown, C. L. Davis 
University of Louisville, Louisville, KY 40292, USA 

J. Allison, N. R. Barlow, R. J. Barlow, P. A. Hart, M. C. Hodgkinson, G. D. Lafferty, A. J. Lyon, 

J. C. WiUiams 

University of Manchester, Manchester M13 9PL, United Kingdom 

A. Farbin, W. D. Hulsbergen, A. Jawahery, D. Kovalskyi, C. K. Lae, V. Lillard, D. A. Roberts 
University of Maryland, College Park, MD 20742, USA 

G. Blaylock, C. Dallapiccola, K. T. Flood, S. S. Hcrtzbach, R. Kofler, V. B. Koptchev, T. B. Moore, 

S. Sarcmi, H. Stacnglc, S. Willocq 

University of Massachusetts, Amherst, MA 01003, USA 
^Now at Department of Physics, University of Warwick, Coventry, United Kingdom 



4 



R. Cowan, G. ScioUa, S. J. Sckula, F. Taylor, R. K. Yamamoto 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, MA 02139, USA 

D. J. J. Mangeol, P. M. Patel, S. H. Robertson 
McGill University, Montreal, QC, Canada H3A 2T8 

A. Lazzaro, V. Lombardo, F. Palombo 
Universitd di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 

J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Krocger, J. Reidy, D. A. Sanders, D. J. Summers, 

H. W. Zhao 

University of Mississippi, University, MS 38611, USA 

S. Brunct, D. Cote, P. Taras 
Universite de Montreal, Laboratoire Rene J. A. Levesque, Montreal, QC, Canada H3C 3J7 

H. Nicholson 

Mount Holyoke College, South Hadley, MA 01075, USA 

N. Cavallo,^ F. Fabozzi,^ C. Gatto, L. Lista, D. Monorchio, P. Paolucci, D. Piccolo, C. Sciacca 
Universitd di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 

M. Baak, H. Bulten, G. Raven, H. L. Snoek, L. Wilden 

NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, 

The Netherlands 

C. P. Jessop, J. M. LoSecco 
University of Notre Dame, Notre Dame, IN 46556, USA 

T. Allmendinger, K. K. Gan, K. Honscheid, D. Hufnagel, H. Kagan, R. Kass, T. PuUiam, A. M. Rahimi, 

R. Ter-Antonyan, Q. K. Wong 

Ohio State University, Columbus, OH 43210, USA 

J. Brau, R. Frey, O. Igonkina, C. T. Potter, N. B. Sinev, D. Strom, E. Torrence 
University of Oregon, Eugene, OR 91403, USA 

F. Colecchia, A. Dorigo, F. Galeazzi, M. Margoni, M. Morandin, M. Posocco, M. Rotondo, F. Simonetto, 

R. Stroili, G. Tiozzo, C. Voci 

Universitd di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

M. Benayoun, H. Briand, J. Chauveau, P. David, Ch. de la Vaissiere, L. Del Buono, O. Hamon, 
M. J. J. John, Ph. Leruste, J. Malcles, J. Ocariz, M. Pivk, L. Rods, S. T'Jampens, G. Therin 

Universites Paris VI et VII, Laboratoire de Physique Nucleaire et de Hautes Energies, F-75252 Paris, 

France 

P. F. Manfredi, V. Re 
Universitd di Pavia, Dipartimento di Elettronica and INFN, 1-27100 Pavia, Italy 
^Also with Universita della Basilicata, Potenza, Italy 



5 



p. K. Bchcra, L. Gladney, Q. H. Guo, J. Panetta 
University of Pennsylvania, Philadelphia, PA 19104, USA 



C. Angelini, G. Batignani, S. Bettarini, M. Bondioli, F. Bucci, G. Calderini, M. Carpinelli, F. Forti, 
M. A. Giorgi, A. Lusiani, G. Marchiori, F. Martinez- Vidal,^ M. Morganti, N. Neri, E. Paoloni, M. Rama, 

G. Rizzo, F. Sandrelli, J. Walsh 

Universita di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 

M. Haire, D. Judd, K. Paick, D. E. Wagoner 
Prairie View A&M University, Prairie View, TX imS, USA 

J. Biesiada, N. Danielson, P. Elmer, Y. P. Lau, C. Lu, V. Miftakov, J. Olsen, A. J. S. Smith, A. V. Telnov 

Princeton University, Princeton, NJ 08544, USA 

F. Bellini, G. Cavoto,^ R. Faccini, F. Ferrarotto, F. Ferroni, M. Gaspero, L. Li Gioi, M. A. Mazzoni, 
S. Morganti, M. Pierini, G. Piredda, F. Safai Tehrani, C. Voena 

Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

S. Christ, G. Wagner, R. Waldi 
Universitdt Rostock, D-18051 Rostock, Germany 

T. Adye, N. De Greet, B. Franek, N. I. Geddes, G. P. Gepal, E. O. Olaiya 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

R. Aleksan, S. Emery, A. Gaidot, S. F. Ganzhur, P.-F. Giraud, G. Hamel de Monchenault, W. Kozanecki, 
M. Legendre, G. W. London, B. Mayer, G. Schott, G. Vasseur, Ch. Yeche, M. Zito 

DSM/Dapnia, CEA/Saclay, F-91191 Gif-sur-Yvette, France 

M. V. Purohit, A. W. Weidemann, J. R. Wilson, F. X. Yumiceva 
University of South Carolina, Columbia, SC 29208, USA 

D. Aston, R. Bartoldus, N. Berger, A. M. Boyarski, O. L. Buchmueller, R. Glaus, M. R. Convery, 
M. Cristinziani, G. Dc Nardo, D. Dong, J. Dorfan, D. Dujmic, W. Dunwoodie, E. E. Elsen, S. Fan, 
R. C. Field, T. Glanzman, S. J. Gowdy, T. Hadig, V. Halyo, C. Hast, T. Hryn'ova, W. R. Innes, 
M. H. Kelsey, P. Kim, M. L. Kocian, D. W. G. S. Leith, J. Libby, S. Luitz, V. Luth, H. L. Lynch, 
H. Marsiske, R. Mcssncr, D. R. MuUcr, C. P. O'Grady, V. E. Ozcan, A. Pcrazzo, M. Perl, S. Pctrak, 
B. N. Ratcliff, A. Roodman, A. A. Salnikov, R. H. Schindler, J. Schwiening, G. Simi, A. Snyder, A. Soha, 
J. Stelzer, D. Su, M. K. Sullivan, J. Va'vra, S. R. Wagner, M. Weaver, A. J. R. Weinstein, 
W. J. Wisniewski, M. Wittgen, D. H. Wright, A. K. Yarritu, C. C. Young 

Stanford Linear Accelerator Center, Stanford, CA 94309, USA 

P. R. Burchat, A. J. Edwards, T. \. Meyer, B. A. Petersen, C. Roat 
Stanford University, Stanford, CA 94305-4060, USA 

S. Ahmed, M. S. Alam, J. A. Ernst, M. A. Saeed, M. Saleem, F. R. Wappler 

State University of New York, Albany, NY 12222, USA 

''Also with IFIC, Institute dc Fi'sica Corpuscular, CSIC-Universidad de Valencia, Valencia, Spain 
^Also with Princeton University, Princeton, USA 



6 



W. Bugg, M. Krishnamurthy, S. M. Spanier 
University of Tennessee, Knoxville, TN 37996, USA 

R. Eckmann, H. Kim, J. L. Ritchie, A. Satpathy, R. F. Schwitters 
University of Texas at Austin, Austin, TX 78712, USA 

J. M. Izen, I. Kitayama, X. C. Lou, S. Ye 
University of Texas at Dallas, Richardson, TX 75083, USA 

F. Bianchi, M. Bona, F. Gallo, D. Gamba 
Universita di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

L. Bosisio, C. Cartaro, F. Cossutti, G. Delia Ricca, S. Dittongo, S. Grancagnolo, L. Lanceri, P. Poropat,^ 

L. Vitale, G. Vuagnin 

Universita di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 

R. S. Panvini 
Vanderbilt University, Nashville, TN 37235, USA 

Sw. Banerjee, C. M. Brown, D. Fortin, P. D. Jackson, R. Kowalewski, J. M. Roney, R. J. Sobie 
University of Victoria, Victoria, BC, Canada V8W 3P6 

H. R. Band, B. Cheng, S. Dasu, M. Datta, A. M. Eichenbaum, M. Graham, J. J. Hollar, J. R. Johnson, 
P. E. Kutter, H. Li, R. Liu, A. Mihalyi, A. K. Mohapatra, Y. Pan, R. Prepost, P. Tan, J. H. von 

Wimmersperg-Toeller, J. Wu, S. L. Wu, Z. Yu 

University of Wisconsin, Madison, WI 53706, USA 

M. G. Greene, H. Neal 
Yale University, New Haven, CT 06511, USA 



® Deceased 



7 



1 Introduction 



In the Standard Model (SM), CP violation arises from a single phase in the three- generation 
Cabibbo-Kobayashi-Maskawa quark-mixing matrix [1]. Possible indications of physics beyond the 
SM may be observed in time-dependent CP asymmetries of B decays dominated by penguin- 
type diagrams to states such as r]'K^, K+R-R^, and /o(980)ir° [2]. Neglecting CKM- 
suppressed amplitudes, these decays carry the same weak phase as the decay B'^ — J/'^R^ [3]. 
As a consequence, their mixing-induced CP-violation parameter is expected to be —r]f x sin2/3 = 
-r]f X 0.74 zb 0.05 [4] in the SM, where /3 = arg [ -T^^T/'^/y^^y^^ ] and r?/ is the CP eigenvalue of 
the final state /, which is for /o(980)if°. There is no direct CP violation expected in these 
decays since they are dominated by a single amplitude in the SM. Due to the large virtual masses 
occurring in the penguin loops, additional diagrams with non-SM heavy particles in the loops and 
new CP-violating phases may contribute. Measurements of CP violation in these channels and 
their comparisons with the SM expectation are therefore sensitive probes for physics beyond the 
SM. 

We present the preliminary results of an update of a measurement [5] of CP-violating asymme- 
tries in the penguin-dominated decay B^ fo^g ^ from a time-dependent maximum-likelihood 
analysis. We restrict the analysis to the region of the 7r+7r~Xg Dalitz plot that is dominated by 
the /o and we refer to this as the quasi-two-body (Q2B) approach. Effects due to the interference 
between the /o and the other resonances in the Dalitz plot are taken as systematic uncertainties. 

The structure of the scalar meson /o has been discussed for decades and is still obscure. There 
were attempts to interpret it as RR molecular states [6], four-quark states [7] and normal qq 
states [8]. However, recent studies of ^ 7/0 (fo 77) [9, 10] and Df Jqtt^ [11] decays 
favor the qq state models. In this interpretation the flavor content of the /q is given by /o = 
cos(^s)ss -f- sin(0s)nn, with nn = {uu + dd)/\/2. A mixing phase of 0s = —48° ± 6° has been 
experimentally determined from — > 7/0 decays [10]. If the assumption is true that the /o state 
has a sizable content of ss, then the decay B^ fo^^s would be dominated by the penguin 
transition, b sss (cf. Fig. 1(b)). Thus, we expect that a measurement of mixing-induced CP 
violation leads to >§ ~ — sin2/?, where S is the coefficient of the sine modulation term [2]. 

The data used in this analysis were accumulated with the BABMl detector [12] at the PEP-II 
asymmetric-energy e+e^ storage ring at SLAC. The data sample consists of an integrated luminosity 
of 192 fb^^ collected at the ^(45") resonance ( "on-rcsonance" ) corresponding to 209 x 10^ BB pairs, 
and 11.8 fb"^ collected about 40MeV below the T{AS) ("off-resonance"). In Ref. [12] we describe 

'''Throughout the paper fo refer to the /o(980) and its decay to tt+tt". In addition, charge conjugate decay modes 
are assumed unless exphcity stated. 




Figure 1: The color- suppressed tree (a) and dominant gluonic penguin (h) are diagrams that could 
contribute to the decay B^ — > /o(980)i(rg. 
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the silicon vertex tracker and drift chamber used for track and vertex reconstruction, and the 
Cherenkov detector (DIRC), the electromagnetic calorimeter (EMC), and the instrumented flux 
return (IFR) used for particle identification. 

If we denote by At the difference between the proper times of the decay of the fully reconstructed 
S° foK^ (^OgJ and the decay of the other meson (-Btag)"' time-dependent decay rate /Qt^g 
is given by 

e-\^t\/r / \ 
/Qtag (At) = ( 1 + QtagS sin(AmdAt) - QtagC cos(AmdAt) 1 (1) 

where Qtag = 1(~1) when the tagging meson -Btag ^ B^(B^), r is the mean lifetime, and 
Am^ is the B^B^ oscillation frequency corresponding to the mass difference. The parameter S is 
non-zero if there is mixing-induced CP violation, while a non-zero value for C would indicate direct 
CP violation. 



2 Candidate Selection 

We reconstruct B^ — *■ /o(— Tr+Tr")^^^ candidates from combinations of two tracks and a de- 
caying to 7r"'"7r~. For the tt'^tt" pair from the /o candidate, we use information from the tracking 
system, EMC, and DIRC to remove tracks consistent with electron, kaon, and proton hypotheses. 
In addition, we require at least one track to have a signature in the IFR that is inconsistent with 
the muon hypothesis. The mass of the /o candidate must satisfy 0.86 < m(7r"'"7r~) < l.lOGeV/c^. 
To reduce combinatorial background from low energy pions, we require | cos 0(7r''")| < 0.9, where 
^(vr"'") is the angle between the positive pion direction in the /o rest frame and the /o flight direc- 
tion in the laboratory frame. The candidate is required to have a mass within lOMeV/c^ of 
the nominal mass [14] and a decay vertex separated from the B^ decay vertex by at least five 
standard deviations. In addition, the cosine of the angle between the flight direction and the 
vector between the /o and the vertices must be greater than 0.99. 

Two kinematic variables are used to discriminate between signal- decays and combinatorial 
background. One variable is the difference A£' between the measured center-of-mass (CM) energy 
of the B candidate and sfsjl^ where -y/s is the CM energy. The other variable is the beam- 
energy-substituted mass mES = \j (s/2 + Pi • tpbY /E'f — p^, where the B momentum ps and 
the four-momentum of the initial state {Ei, pj) are defined in the laboratory frame. We require 
5.23 < mES < 5.29GeV/c2 and \/\E\ < 0.1 GeV. 

Continuum e^e~ qq{q = u, d, s, c) events are the dominant background. To enhance discrim- 
ination between signal and continuum, we use a neural network (NN) to combine four variables: 
the cosine of the angle between the B^^^ direction and the beam axis in the CM, the cosine of the 
angle between the thrust axis of the B*^^^ candidate and the beam axis, and the zeroth and second 
angular moments Lo,2 of the energy flow about the iJ^ec thrust axis. The moments are defined by 
Lj = Y^- Pi X I cos 9i p , where 6i is the angle with respect to the B^^ thrust axis of the track or 
neutral cluster i, and pi is its momentum. The sum excludes the B^^ candidate. The thrust axis is 
defined as the direction that maximizes the sum of the longitudinal momenta of the B^^^ daughters. 
The NN is trained with off-resonance data and simulated signal events. The final sample of signal 
candidates is selected with a cut on the NN output > —1.5, which retains ~ 97% and 52% of the 
signal and continuum, respectively. 

"The B°ag is so called because its flavor is determined using the tagging algorithm of Ref. [13]. 
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The signal efficiency determined from Monte Carlo (MC) simulation is (38.7 ± 0.4)%. MC 
simulation shows that 4.7% of the selected signal events are mis-reconstructed, mostly due to 
combinatorial background from low-momentum tracks used to form the /o candidate. In total, 
12586 on-resonance data events pass all selection criteria. 

3 Background from other B Decays 

We use MC-simulated events to study the background from other B decays. The charmless decay 
modes are grouped into eight classes with similar kinematic and topological properties. The modes 
that decay to the Tr~^Tr~K^ final state are of particular importance since they have signal-like AE 
and rriES distributions and their decay amplitudes interfere with the foK^ decay amplitude. Among 
these modes are p'^K^, fo{1370)Kg, /2(1270)ifg, K*+7r'~ (including other kaon resonances decaying 
to -K'°7r+), and non-resonant 7t~^7t Kg decays. The mode p^K^ is particularly important because it 
has r]f = —1 and thus any p^K^ events misidentified as signal will dilute the observed CP asymmetry 
in our data. The inclusive charmless iT'^iT~Kg branching fraction (23.4 ± 3.3) x 10~^ [4], together 
with the available exclusive measurements [4], are used to infer upper limits on the branching 
fractions of these decays. Along with selection efficiencies obtained from MC, these branching 
fractions are used to estimate the expected background. The charmed decays — > D~Tr~^ — > 
KgTr~7r~^ and 5+ —>■ iD°7r+ — K^tt^tt''" contribute significantly to the selected data sample. Each 
of these modes is treated as a separate class. Two additional classes account for the remaining 
neutral and charged 6 — c decays. In the selected data sample we expect 106 it 23 charmless and 
218 ± 93 6 ^ c events. 



4 Meiximum-Likelihood Fit 

The time diff'erence At is obtained from the measured distance between the z positions (along 
the beam direction) of the iJ^ec -^tag decay vertices, and the boost fi^ = 0.56 of the e+e" 
system [13, 15]. To determine the flavor of the B^^^^ we use the tagging algorithm of Ref. [13]. 
This produces four mutually exclusive tagging categories. We also retain untagged events in a fifth 
category to improve the efficiency of the signal selection. 

We use an unbinned extended maximum-likelihood fit to extract the /o-K's event yield, the CP 
parameters defined in Eq. (1), and the /o resonance parameters. The likelihood function for the 
Nk candidates tagged in category k is 



Ck = e-^^n ^sek [(1 - /W)^^ + /^R^^-^^] + NckV^k + E (2) 

i=l ^ j=l ^ 



where Nj^ is the sum of the signal, continuum and the ns -B-background yields tagged in category 
k, Ng is the number of foKg signal events in the sample, is the fraction of signal events tagged in 
category k, f^-^ is the fraction of mis-reconstructed signal events in tagging category k, Nc^k is the 
number of continuum background events that are tagged in category A;, and NB,j^j,k is the number 
of i?-background events of class j (sec section 3) that are tagged in category k. The i3-background 
event yields arc fixed parameters, with the exception of the D~tt~^ yield. Since B^ D^tt^ events 
have a characteristic distribution in cos0(7r+), well separated from continuum and foKg events, 
the D~Tr~^ is free to vary in the fit along with the signal and continuum yields. The total likelihood 
£. is the product of the likelihoods for each tagging category. 
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The probability density functions (PDFs) V^~^^, V^~^^, and Vf,^, for correctly recon- 
structed signal, mis-reconstructed signal, continuum background and B-background class j, respec- 
tively, are the products of the PDFs of six discriminating variables. The correctly reconstructed sig- 
nal PDF is thus given by: "^^ = P^-^^{mEs) - V^-^^iAE) - V^'^^iNN) ■V'^-^^{\ cos6'(7r+)|) • 
V^~^^{m{7r+TT-)) ■ Pf"'^^(At), where V^''~^^{At) contains the time-dependent CP parameters 
defined in Eq. (1), diluted by the effects of mis-tagging and the At resolution. 

The fractions of mis-reconstructed signal events in each tagging category are estimated by MC 
simulation. The uies, AE, NN, | cos 0(7r"'")|, and m(7r"'"7r~) PDFs for signal and B background 
are taken from the simulation except for the means of the signal Gaussian PDFs for ttt-es and AE 
as well as the mass and width of the /o, which are free to vary in the fit. We use a relativistic 
Breit-Wigner function to parameterize the /o resonance. The At-resolution function for signal and 
B-background events is a sum of three Gaussian distributions, with parameters determined by a fit 
to fully reconstructed decays [13]. The continuum At distribution is parameterized as the sum 
of three Gaussian distributions with two distinct means and three distinct widths, which are scaled 
by the At per-event error. For the S-background modes that are CP eigenstates, the parameters C 
and S are fixed to and ib sin 2/3, respectively, depending on their CP eigenvalues. For continuum, 
four tag asymmetries and the five yields Nc\k are free. The signal yield, S, C, and the /o mass 
and width are among the 41 parameters that are free to vary in the fit. The majority of the free 
parameters are used to describe the shape of the continuum background. 

5 Systematic Errors 

The contributions to the systematic error on the signal parameters are summarized in Table 1. To 
estimate the errors due to the fit procedure, we perform fits on a large number of MC samples with 
the proportions of signal, continuum and 5-background events measured from data. Biases of a 
few percent observed in these fits are due to imperfections in the likelihood model such as neglected 
correlations between the discriminating variables of the signal and i?-background PDFs and are 
assigned as a systematic uncertainty of the fit procedure. The error due to the fit procedure includes 
these biases added in quadrature with their statistical errors. The expected event yields from the B- 
background modes are varied according to the uncertainties in the measured or estimated branching 
fractions. Since -B-background modes may exhibit CP violation, the corresponding CP parameters 
are varied within their physical ranges. We vary the parameters of the At model and tagging 
fractions incoherently within their errors and assign the observed changes, summed in quadrature, 
as a systematic error. The uncertainties due to the simulated signal PDFs are obtained from a 



Table 1: Summary of systematic uncertainties. 



Error Source 


S 


c 


Fitting Procedure 


0.06 


0.10 


B-background 


0.04 


0.08 


At Model 


0.01 


0.01 


Tagging 


0.02 


0.01 


Signal Model 


0.02 


0.02 


DCS Decays 


0.01 


0.04 


Am,d and r 


0.00 


0.01 


Q2B Approximation 


0.04 


0.07 


Sub-total 


0.10 


0.15 
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control sample of fully reconstructed D~{-^ K^n~)Tr~^ decays. The systematic errors due to 

interference between the doubly-Cabibbo-suppressed (DCS) b ucd amplitude with the Cabibbo- 
favored 6 cud amplitude for tag-side B decays have been estimated from simulation by varying 
freely all relevant strong phases [16]. The errors associated with Amd and r are estimated by 
varying these parameters within the errors on the world average [14]. 

The systematic error introduced in the Q2B approximation by ignoring interference effects 
between the /o and the other resonances in the Dalitz plot is estimated from simulation by varying 
freely all relative strong phases and taking the largest observed change in each parameter as the 
error. Eleven resonances are used in this study including the three lowest lying p resonances, 
/o(980), /o(1370), /2(1270), and the K*^ and higher kaon states. In addition, a non-resonant 
component is allowed. The proportion of each contribution is estimated using known exclusive 
measurements and the inclusive 7r"'"7r~i^^^ rate. The systematic effects due to interference are small 
compared with the statistical error for S and C. 

6 Fit Results 

The maximum-likelihood fit results in the CP-violation parameters: 

S = -0.95ig:|±0.10, 
C = -0.24 ±0.31 ±0.15, 

where the first errors are statistical and the second are systematic. The improvement in the error 
with respect to the previous result (127 x 10^ T{4S) ^ BB decays, (Tstat(5') =-o'li) mainly 
to the increased luminosity, but is due also to one event with large signal probability and to the 
proximity of the measured S and the physical limit (IS*] < 1). We find an foKg event yield of 
152.4 ± 18.5 which is consistent with the previously measured branching fraction [5]. 

Figure 2 shows distributions of AE, niES, \ cos ^(7r+)| and m(7r"'"7r~), that are enhanced in signal 
content by cuts on the signal-to-continuum likelihood ratios of the other discriminating variables. 
The time-dependent distributions and asymmetry A^o/^o ~ (-^3° ~ ^B°)/i-^B° + ^b") 
tagged events are presented in Fig. 3. 

We validated the stability of the nominal fit by testing different fit configurations where each 
configuration had a discriminating variable removed. As another cross-check, we allow the B^ 
lifetime, r^o, to vary. We find r^o = (1.52 ± 0.22) ps, in agreement with the world average r^o = 
(1.536 ± 0.014) ps [4], and the remaining free parameters are consistent with the nominal fit. 

7 Summary 

In summary, we have presented an updated preliminary measurement of the CP-violating asym- 
metries in B^ /o(980)(^ Tr+Tr")^^ decays. Our results for S and C are consistent with the 
Standard Model. The hypothesis of no mixing-induced CP violation is excluded at the 2.3cr level. 
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Figure 2: Distributions of (clockwise from top left) AE, rriEs, |cos0(7r+)|, m(7r"'"7r~) and the 
NN output for samples enhanced in foK^ signal (purity is ~ 45%.) The solid curve represents 
a projection of the maximum-likelihood fit result. The dashed curve represents the contribution 
from continuum events, and the dotted line (middle) indicates the combined contributions from 
continuum events and B backgrounds. For presentation purposes, the region 0.765 < | cos0(7r+)| < 
0.81 has been removed to suppress the contribution from D~Tr~^ events. 
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Figure 3: The signal enhanced time distributions tagged as B^^^ (top) and B^^^ (middle), and the 
asymmetry, A^o/bo (bottom). The solid curve is a projection of the fit result. The dashed line 
is the distribution for continuum background and the dotted line is the total B- and continuum- 
background contribution. 



14 



8 Acknowledgments 



We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have made this work possible. The success of 
this project also relies critically on the expertise and dedication of the computing organizations 
that support BABMl. The collaborating institutions wish to thank SLAC for its support and the 
kind hospitality extended to them. This work is supported by the US Department of Energy and 
National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), 
Institute of High Energy Physics (China), the Commissariat a I'Energie Atomique and Institut 
National de Physique Nucleaire et de Physique des Particules (France) , the Bundesministerium fiir 
Bildung und Forschung and Deutsche Forschungsgcmcinschaft (Germany), the Istituto Nazionale 
di Fisica Nuclearc (Italy), the Foundation for Fundamental Research on Matter (The Netherlands), 
the Research Council of Norway, the Ministry of Science and Technology of the Russian Federation, 
and the Particle Physics and Astronomy Research Council (United Kingdom). Individuals have 
received support from CONACyT (Mexico), the A. P. Sloan Foundation, the Research Corporation, 
and the Alexander von Humboldt Foundation. 

References 

[1] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. Theor. 
Phys. 49, 652 (1973). 

[2] Y. Grossman and M. P. Worah, Phys. Lett. B 395, 241 (1997); M. Ciuchini et ai, Phys. Rev. 
Lett. 79, 978 (1997); D. London and A. Soni, Phys. Lett. B 407, 61 (1997). 

[3] R. Fleischer, Z. Phys. C 62, 81 (1994); Y. Grossman, Z. Ligeti, Y. Nir and H. Quinn, Phys. 
Rev. D 68, 015004 (2003). 

[4] The Heavy Flavor Averaging Group (HFAG), J. Alexander et ai, 
http:/ /www.slac.stanford.edu/xorg/hfag/ 

[5] BABAR Collaboration, B. Aubert et al, hep-ex/0406040 (2004). 

[6] J. Weinstein and N. Isgur, Phys. Rev. Lett. 48, 659 (1982); Phys. Rev. D27, 583 (1983); Phys. 
Rev. D41, 2236 (1990); M.P. Locher et al, Eur. Phys. J. C4 317 (1998). 

[7] R.J. Jaffe, Phys. Rev. D15, 267 (1997); M. Alford and R.J. Jaffe, Nucl. Phys. B578, 367 
(2000). 

[8] N.A. Torinqvist, Phys. Rev. Lett. 49, 624 (1982); N.A. Tornquits and M. Roos, Phys. Rev. 
Lett. 76, 1575 (1996). 

[9] F. De Fazio and M.R. Pennington, Phys. Lett. B521, 15 (2001); R. Delborgo, D. Liu and 
M.D. Scadron, Phys. Lett. B446, 332 (1999); T.M. Aliev et al, Phys. Lett. B527, 193 (2002). 

[10] A.V. Anisovich, V.V. Anisovich and V.A. Nikonov, hep-ph/0011191 (2000) 

[11] F. Kleefeld et al, Phys. Rev. D66, 034007 (2002); E. van Beveren, G. Rupp and M.D. Scadron, 
Phys. Lett. B495, 300 (2000); A. Deandrea et al, Phys. Lett. B502, 79 (2001). 



15 



[12] BABAR Collaboration, B. Aubert et al, Nucl. Instrum. Methods Phys. Res., Sect. A 479, 1 
(2002). 

[13] BABAR Collaboration, B. Aubert et al, Phys. Rev. D66, 032003 (2002). 

[14] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 

[15] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 89, 281802 (2002). 

[16] O. Long and M. Baak, R.N. Cahn, D. Kirkby, Phys. Rev. D 68, 034010 (2003). 

[17] E791 Collaboration, E.M. Aitala et al, Phys. Rev. Lett. 86, 765 (2001). 



16 



